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ABSTRACT. Modeling of the velocity field of the Local Supercluster leads to 
the conclusion that there must be a substantial differential between the mass- 
to-light ratios of most galaxies and those in a small number of special places. 
Specifically, M/L = IOOOA/q/Lq is indicated for the full Virgo Cluster while 
M/L = 150Mq / Lq is appropriate for most environments. It is argued that 
a higher M/L value is characteristic of E/SO knots, regions where the galaxy 
density is high and collapse time scales are short. Regions larger than galaxies 
where collapse is ocurring only on the order of the age of the universe, the 
environment where spirals and irregulars predominate, are characterized by 
a low M/L. Overall, the mean density of the universe is well below closure. 
There could be comparable total mass in the relatively rare E/SO knots as 
there is in the environments dominated by spirals and irregulars. 



1 Introduction 

The concept of 'biasing' arose out 
of a discussion by Kaiser (1984) that 
different classes of astronomical ob- 
jects may trace the underlying mass 
distribution to greater or lesser de- 
grees. It was appreciated, especially 
as comparisons were made with N- 
body simulations, that some fraction 
of mass could be only poorly corre- 
lated with the observable galaxies. A 
simple, perhaps simplistic, descrip- 
tion is provided by the bias parame- 
ter of linear theory, b = 8 g / 5 m , where 
the observed galaxy fluctuation field 
is 8 g and the mass density fluctu- 
ation field is 8 m . The expectation 
has arisen from comparisons between 
models and the observed spectrum 
of irregularities (Davis et al. 1985) 
that mass would be more widely dis- 
tributed than clusters and even field 
galaxies, so b > 1. In this case, it 
would be said that galaxies have a bi- 
ased distribution compared with the 
mass. 

While this sense of biasing may 
well occur, the focus of this paper is 
on the opposite possibility that mat- 
ter may be more concentrated than 
light in some cosmologically impor- 
tant situations. At least locally, it 
could be that b < 1, a condition we 
will call 'antibiasing'. 



2 The Mean Density of the 
Universe 

Dynamical studies of galaxy flows 
can give estimates of the mean den- 
sity (Dekel 1994; Strauss & Willick 
1995). There is a dichotomy of re- 
sults. One method compares the di- 
vergence of the velocity field, which 
gives the mass density field, with 
some catalog of the observed dis- 
tribution of galaxies. This method, 
called POTENT, gives values of (3 = 
S1q' 6 /6 that are indicative of a uni- 
verse near closure density: f3j = 
0.89 ± 0.24 (Sigad et al. 1998) and 
p = 0.74±0.26 (Hudson et al. 1995) 
where the subscripts / and O indi- 
cate comparison with an Infrared As- 
tronomy Satellite (IRAS) and opti- 
cal catalog, respectively. Baker et al. 
(1998) find 6 /6/ ~ 1.4, so these re- 
sults would be compatible with flo ~ 
1 and bj ~ 1. The uncertainties are 
95% confidence levels. 

On the other hand, methods that 
compare expectation and observed 
velocities find lower values of f3. Here, 
both linear and non-linear studies 
reach similar conclusions. The lin- 
ear approximation is valid if 8 m < 1 
(Peebles 1980): 

v=£ / d 3 r A(Q(r'-r) 
4n J | r' - r | 3 . 
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(1992) found O < 0.25 (ft < 
0.35), Reiss et al. (1997) found ft = 

0A - a25 ' da Costa et al - (!998) found 
ft = 0.6 ±0.2, and Willick & Strauss 
(1998) found ft = 0.50±0.10. A non- 
linear methodology we call 'Least 
Action' has been described by Pee- 
bles (1989, 1990, 1994, 1995) and 
Shaya, Peebles & Tully (1995). The 
non-linear studies determine, not (3, 
but constraints on the parameter 
space (M/L, to), where M/L is the 
conversion ratio from observed lumi- 
nosities to mass and to is the age of 
the universe. The M/L values that 
are determined can be normalized 
by the mean luminosity density of 
the universe to constrain something 
equivalent to /3. The results from our 
most recent analysis (Tully 1997) can 
be translated as fi = 0.31 ± 0.14 
(ft = 0.43 ± 0.20). All the linear 
and non-linear studies that compare 
model velocities to observed veloci- 
ties are in reasonable agreement and 
indicate that if Qq = 1 then there 
needs to be considerable bias, b; ~ 2. 



3 Least Action 

Shaya ct al. (1995) describe our 
methodology for the reconstruction 
of the orbits of galaxies and groups 
using a variational principle. Tully 
(1997) gives updated results. Galax- 
ies follow paths that extremize the 
'action', the integral of the La- 
grangian through time. In its sim- 
plest incarnation, the model can be 
restricted to only two free parame- 
ters. We have a catalog of galaxies 
within a velocity limit (3000 km s _1 ) 
with known positions on the sky, red- 
shifts, and luminosities, so one of our 
free parameters, the relation between 
mass and light (M/L), specifics the 
mass associated with galaxies where 
they are today in redshift space. The 
other free parameter is the age of the 
universe, to, that specifies how long 
these masses have been tugging on 
each other. At the next level of com- 
plexity, a third parameter is added 
to take some account of biasing. A 
softening parameter is introduced to 
modify the force law so forces are re- 
duced at close range to mimic the ef- 
fects of inter-penetrating halos. 

Once the three free parameters 
have been fixed, a Least Action re- 
enactment of the constituent orbits 



gets all the galaxies to the right 
positions on the sky and the right 
rcdshifts (as boundary conditions) 
and specifies the other three ele- 
ments of phase-space. The specifi- 
cation of model distances is what 
we care about because measured dis- 
tances exist for ~ 1/3 of the ob- 
jects. A x 2 cvaluator is formulated 
to compare model and observed dis- 
tance moduli, fi mt i and fi 0} i respec- 
tively, for the i th object: 



Uog(cz 0>i /cz m>i ) s 



+( 



Figure 1 illustrates the \ 2 val- 
ues in the domain of the two princi- 
pal free parameters (M/L, to) with 
the softening parameter fixed at 
300 km s _1 . (Note: since this model 
was created we have revised our ta- 
ble of distances; the Hubble Space 
Telescope observations of cepheids in 
external galaxies has resulted in a 
10% upward revision of distances and 
ages: see Tully 1998). It can be seen 
that models with Ho substantially 
less than unity are favored. Mod- 
els with M/L ~ 200M Q /L Q , t ~ 
10 Gyr conform to observations. 



4 The Virgo Cluster 

Although the model described 
above provides a good description of 
the velocity field over most of the Lo- 
cal Supcrclustcr, it fails miserably in 
one important region. A substantial 
infall pattern is observed around the 
Virgo Cluster but such a low value 
as M/L ~ 200Mq/Lq cannot ex- 
plain the motions that are seen. It 
hardly helps to simply increase M / L 
for all objects. Doing so effectively 
increases the model Qo, hence re- 
duces to for a given distance scal- 
ing. There is then not enough time 
to build up the observed velocities. 
High global M/L cannot explain the 
infall velocities and the \ 2 figure of 
merit blows up. 

The very strong conclusion is that 
we need a large differential in M/L 
between Virgo and most of the rest 
of the Local Supercluster. We need 
a modest M/L outside Virgo to get 
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Figure 1. Contours of \ 2 as a function of the two free parameters M/L 
and to with the force softening parameter set at 300 km s~ 1 . Contours at lo 
intervals with the dark contour at the 2o level. This model is calibrated with 
distances that do not reflect a 10% increase in the zero-point resulting from 
recent cepheid distance determinations. 



a good x 2 fit overall but we need a 
large M/L in the Virgo Cluster to 
explain the observed infall pattern. 

Figure 2 illustrates the nature of 
our observational constraints. The 
group 11-4 (Tully 1987) is an en- 
tity in the Southern Extension of 
the Virgo Cluster, at an angle of 
9.3° from the cluster center. The 
curves in each panel show the run of 
velocities with distance anticipated 
by two alternative models: on top, 
a model with Q = 0.3, M/L = 
2OOM0 /Lq and, on the bottom, a 
model with H = 0.3, M/L = 
15OM0 /Lq for most galaxies but 
M/L = 1OOOM /L for Virgo and 
a small number of other groups dom- 
inated by ellipticals and SOs. The 
points with error bars locate galax- 
ies in group 11-4 with distance es- 
timates. Clearly, the model that de- 
fines the curve in the top panel can- 
not explain the velocities that are ob- 
served. The model associated with 
the bottom panel does an adequate 
job. 

All galaxies within the infall re- 
gion can be considered in the same 
light. In each case, the run of veloci- 
ties with distance along the pertinent 
line-of-sight can be calculated for any 
model. For a given to, there will be 
a threshold {M / L) minimum that al- 
lows the observed velocity within the 
infall envelope. 

It is helpful to appreciate that we 
do not need very precise distance es- 
timates for this analysis. It is suffi- 



cient to discriminate between three 
possibilities: is a galaxy to the fore- 
ground, within the infall region, or 
in the background? Since the Virgo 
multi-value region is seen to extend 
to ~ 25° from the Virgo core, the in- 
fall region extends to ~ 1 . 2 m in front 
of the mean cluster distance and to ~ 
0.8 m behind the mean distance. Our 
distance estimators are adequate to 
discriminate between the three alter- 
natives. The case of Fig. 2 provides 
an illustration. 

The Virgo Cluster caustic lies at ~ 
6° from the core and the infall is pro- 
nounced out to ~ 20°. Within this 
annulus, we currently have 35 galax- 
ies with distance determinations that 
unambigously put them in the in- 
fall zone. These galaxies are falling 
in a quasi-laminar flow into the clus- 
ter. Figure 3 shows the constraints 
these galaxies provide for the same 
two models considered in Fig. 2. The 
circles show the observed velocities 
and angular distances from the clus- 
ter center of galaxies that are certi- 
fied to be within the Virgo infall zone 
on the basis of distance measure- 
ments. The small open circles cor- 
respond to galaxies projected onto 
the cluster and, hence, are probable 
cluster members. The 35 filled circles 
correspond to galaxies projected out- 
side the cluster proper. The horizon- 
tal bars in each panel indicate the 
maximum and minimum projected 
velocities along lines of sight to in- 
dividual galaxies or groups of galax- 
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Figure 2. Example of model velocities along a line-of-sight close to the di- 
rection of the Virgo Cluster. In this case, the sight line is through Group 
11-4 ■ The points with errors correspond to galaxies in this group with dis- 
tance determinations. The vertical dashed lines bracket the distance of the 
Virgo Cluster, centered in distance and velocity at the large cross. Top panel: 
M / L = 200 for all entities. Bottom panel: M / L = 1000 for Virgo and other 
E/SO knots, otherwise M/L = 150. The second wave in the velocity curve 
beyond Virgo occurs because the line-of-sight passes near another E/SO knot, 
the Virgo W Cluster, Group 11-24- 



ies within the infall region. The 11- 
4 group used as an example above 
is at an angle of 9.3° and the maxi- 
mum and minimum velocities antic- 
ipated by the single value M/L = 
2OOM0/L0 model (top panels Figs. 



2 and 3) are 1037 and 760 km s" 1 , 
while the maximum and minimum 
velocities anticipated by the model 
with M/L = 1OOOM /L e for Virgo, 
M/L = 150Mq/L q for most other 
objects (bottom panels Figs. 2 and 
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Figure 3. Virgo infall constraints from two models. Horizontal bars are line- 
of-sight extrema velocities for a good overall M/L (top) and for a model with 
a lot of extra mass in Virgo (bottom). 



3) arc 1690 and 283 km s _1 . The 
model illustrated in the top panels 
totally fails to predict the observed 
infall pattern while the model dis- 
played in the bottom panels is ad- 
equate. In this case, the Virgo Clus- 
ter is assigned 1.3 X 1O 15 M . This 
mass is higher than the Virial theo- 
rem mass of 0.7 X 10 15 M Q (Tully & 
Shaya 1984), but the infall mass per- 
tains to a scale ~ 3x larger than the 
Virial radius. 

Fifteen years ago, Tully & Shaya 



(1984) had already found this ba- 



sic result that (M/L) v 



has to be 



much higher than (M/L)fi e i^ in or- 
der to get a sensible infall descrip- 
tion. That previous study was based 
on a fully non-linear but spheri- 
cally symmetric model and one could 
wonder if the obvious departure of 
the observed galaxies from a spher- 
ical distribution in the Local Su- 
pcrcluster was affecting the solu- 
tion. The present model provides a 
non-paramctric description of the su- 
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percluster asymmetry (mass is dis- 
tributed like the observed galaxies) 
and returns essentially the same an- 
swer. The infalling galaxies in Fig. 3 
tend to have higher velocities than 
the cluster mean, partially a hap- 
penstance of uneven filling of the 
volume and partly due to angular 
momentum in the infalling popula- 
tion, naturally provided by the Least 
Action reconstruction. This veloc- 
ity asymmetry was seen in the 1984 
study but could not be described by 
the spherically symmetric model (see 
also Tully & Shaya 1994; Tully 1997). 

5 High M/L in a Small E/SO 
Knot 

Throughout the Local Superclus- 
ter there are places where the galaxy 
density is as high as in the core of 
the Virgo Cluster, though the to- 
tal number of galaxies involved can 
be small. Usually in these cases the 
majority of the galaxies arc ellip- 
ticals and SOs. There is some evi- 
dence that M/L values can be very 
high in these places. A dramatic 
case is the group around NGC 1407, 
Group 51-8 in our catalog (Tully 
1987). Its unusual properties have 
been discussed by Gould (1993) and 
Quintana, Fouque, & Way (1994). 
There are 19 galaxies with known 
rcdshifts associated with the group. 
Only NGC 1407 is an L* galaxy. 
The second most luminous galaxy, 
NGC 1400, is more than three times 
fainter. All but one of the dozen that 
have been typed are E/SO/Sa. Quin- 
tana et al. find the velocity disper- 
sion to be 416^gg kms -1 , which 
gives (M/L) virial = 6OOM /L af- 
ter a tiny correction to our distance 
scale. 

The SO galaxy NGC 1400 de- 
serves special attention. This galaxy 
is blueshifted by 1066 km s _1 with 
respect to the group mean and would 
carry the vast majority of the ki- 
netic energy of the group (~ 85% of 
the kinetic energy follows from the 
line-of-sight velocities). One way to 
rationalize this unusual situation is 
to suppose that NGC 1400 and the 
other galaxies carry only a small frac- 
tion of the mass of the gravitational 
potential well. 



6 Other E/SO Knots in the Lo- 
cal Supercluster 

In a complete volume-limited sam- 
ple in our vicinity of the Local Su- 
percluster (extending to 25h^ Mpc 
with distances given by a model with 
spherical Virgo infall), 69% of the 
galaxies, and 77% of the light in 
galaxies, are in 179 groups (Tully 
1987). Spiral and irregular galax- 
ies predominate in the overwhelm- 
ing majority of these groups. In these 
most common cases, the E/SO rep- 
resentation is inevitably less than 
20%, 1-dimensional group velocity 
dispersions arc ~ 100 km s , and 
group crossing-times arc ~ 40% of 
the Hubble time. However, eight 
of these groups are quite distinct. 
Indeed, they arc among the best 
known nearby groups because of 
their compact nature. In these eight 
cases, more than half of the galax- 
ies are E/SO. These early type galax- 
ies are inevitably restricted to vol- 
umes only a few hundred kpc across. 
The groups usually have higher ve- 
locity dispersions than their spiral- 
dominated counterparts and crossing 
times are < 10% of the Hubble time. 
The distinction between spiral-rich 
and E/SO-rich groups is quite pro- 
nounced, obscured only partially by 
the not-infrequent incidence of spi- 
rals at the periphery of E/SO knots. 

Table 1 provides a list of the com- 
plete sample of nearby groups dom- 
inated by elliptical and SO galax- 
ies (the asterisk in the final col- 
umn indicates inclusion in the com- 
plete sample), plus a few other such 
groups at greater distances but still 
within 3000 kms -1 . The most lu- 
minous and the faintest groups on 
the list have already been discussed. 
The Virgo Cluster and the NGC 1407 
Group span a range from two dozen 
L* galaxies and 10 15 Mq to a sin- 
gle L* galaxy and 2 X 10 13 M Q . For- 
nax, Coma I, Antlia, and NGC 1566 
all have internal velocity dispersions 
that hint at large M/L values (Tully 
1987). 

7 Summary 

The evidence for large M/L varia- 
tions comes from combining two dis- 
tinct sources of dynamical informa- 
tion. On the one hand, the Least Ac- 
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Tabic 1. E/SO Knots within 3000 km s~ 



Name 


Cat: 


ilog 


V LG 


Log L B 





Virgo 


11 


-1 


986 


12.14 


* 


Fornax 


51 


-1 


1466 


11.52 


* 


Virgo W 


11- 


24 


2225 


11.43 




Coma I 


14 


-1 


983 


11.34 


* 


NGC 5044 


11- 


31 


2559 


11.19 




Antlia 


31 


-2 


2519 


11.12 




Eridanus 


51 


-4 


1415 


11.06 


* 


NGC 1566 


53 


-1 


1076 


11.03 


* 


NGC 4125 


12 


-5 


1436 


11.00 


* 


NGC 5846 


41 


-1 


1810 


10.99 




Leo 


15 


-1 


722 


10.73 


* 


NGC 1407 


51 


-8 


1522 


10.58 


* 



tion reconstruction of orbits is telling 
us that most of the galaxies (ie, the 
spirals and irregulars) live in an en- 
vironment with a modest mean den- 
sity. On the other hand, the large 
infall motions into the Virgo Clus- 
ter and the large internal dispersions 
in other E/S0 knots arc indications 
of a lot of mass in these special lo- 
cations. The Least Action dynami- 
cal measurements encompass larger 
scales than virial measurements of 
internal group motions which prob- 
ably explains why the Least Action 
mass estimates are larger by up to a 
factor 2. It is suggested that M/L ~ 
150Mq /Lq for most galaxy associ- 
ations but M/L ~ 1OOOM /L in 
the knots of ellipticals and SOs. Weak 
lensing provides a probe of the dis- 
tribution of matter on superclustcr 
scales for distant objects and it is in- 
teresting that hints are arising from 
that work that mass is strongly 'an- 
tibiased' toward the clusters (Kaiser 
ct al. 1998). 

If E/S0 knots do have much higher 
M/L values than galaxies in lower 
density regions there are several 
potential astrophysical explanations: 
(i) dimming at Lg is anticipated for 
older stellar populations, (ii) strip- 
ping from encounters might shred a 
large fraction of stars into intraclus- 
ter space; Mendez et al. (1997) and 
Ferguson et al. (1998) find enough in- 
tracluster planetary nebulae and red 
giants between the galaxies in Virgo 
to indicate that there are compara- 
ble numbers of stars outside galax- 
ies as in, (iii) dwarf spheroidal galax- 
ies are found in large numbers in 



Virgo (Phillipps ct al. 1998) and For- 
nax (Ferguson & Sandage 1988) but 
evidently not in large numbers in 
Ursa Major or the field (Trentham, 
Tully, & Verheijcn, in preparation) 
(iv) field galaxies may continue to ac- 
crete gas even until today, gas that 
gets converted into stars, while gas 
that arrives in clusters after the ini- 
tial collapse of the cluster is outside 
the Roche-limit of individual galax- 
ies and becomes virializcd as hot gas 
in the intracluster medium (Shaya & 
Tully 1984). It is easy to entertain 
that most or all these effects occur 
and cummulativcly could account for 
an M/L differential of the observed 
factor of 6-7. 

Roughly 8% of the B light within 
3000 km s -1 comes from the galax- 
ies associated with E/S0 knots. How- 
ever, with the M/L differential that 
has been suggested, the E/S0 knots 
could contain ~ 40% of the mass 
directly associated with galaxies. If 
so, there is a dramatic antibiasing of 
matter compared with IRAS-selected 
galaxies or even optically selected 
galaxies. Dark matter is more con- 
centrated in a few places but enough 
so that it would have a cosmological 
impact. In addition, there may also 
be the bias of some matter poorly 
correlated with galaxies. The Least 
Action modeling strongly suggests 
that the overall density of matter is 
modest, U Q ~ 0.25. The POTENT 
observations of f2o = (£>/3) 1,7 can be 
reconciled with a value of Qq well less 
than unity if, overall, bj < 1. 
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